The Edge of Quantum Chaos 
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We identify a border between regular and chaotic quantum 
dynamics. The border is characterized by a power law de- 
crease in the overlap between a state evolved under chaotic 
dynamics and the same state evolved under a slightly per- 
turbed dynamics. For example, the overlap decay for the 
quantum kicked top is well fitted with [l-|-(g-l)(t/r)^]i/'^"«' 
(with the nonextensive entropic index q and r depending on 
perturbation strength) in the region preceding the emergence 
of quantum interference effects. This region corresponds to 
the edge of chaos for the classical map from which the quan- 
tum chaotic dynamics is derived. 

PACS numbers 03.65. Ta, 03.67.-a, 05.20.-y, 05.45. +b 

Classical chaotic dynamics is characterized by strong- 
sensitivity to initial conditions. T-wo initially close points 
move apart exponentially rapidly as the chaotic dynamics 
evolve. The rate of divergence is quantified by the Lya- 
punov exponent [1]. At the border bet-ween chaotic and 
non-chaotic regions (the 'edge of chaos'), the Lyapunov 
exponent goes to zero. Ho-wever, it may be replaced by 
a generalized Lyapunov coefficient [2] describing po-wer- 
la-w, rather than exponential, divergence of classical tra- 
jectories. 

This paper identifies a characteristic signature for the 
edge of quantum chaos. Quantum states maintain a con- 
stant overlap fidelity, or distance, under all quantum dy- 
namics, regular and chaotic. One -way to characterize 
quantum chaos is to compare the evolution of an ini- 
tially chosen state under the chaotic dynamics -with the 
same state evolved under a perturbed dynamics [3] [4] [5] . 
When the initial state is in a regular region of a mixed 
system, a system -with regular and chaotic regions, the 
overlap remains close to one. When the initial state is 
in a chaotic zone, the overlap decay is exponential. This 
paper sho-ws that at the edge of quantum chaos there is 
a region of polynomial overlap decay. 

The Lyapunov exponent description of chaos is as fol- 
lo-ws [1]. If Axq is the distance bet-ween t-wo initial con- 
ditions, -we define ^ — lirriAxo^oi^^) to describe ho-w 
far apart t-wo initially arbitrarily close points become at 
time t. Generally, ^(t) is the solution to the differential 
equation — \i^{t), such that ^(t) = e'^^* (Ai is the 
Lyapunov exponent). When the Lyapunov exponent is 



positive the dynamics described by ^(t) is strongly sensi- 
tive to initial conditions and we have chaotic dynamics. 

This description of chaos -works -well for classical Ne-w- 
tonian mechanics but it cannot hold true for quan- 
tum mechanical -wave functions governed by the linear 
Schrodinger equation. Indeed, like the overlap bet-ween 
Liouville probability densities, the overlap bet-ween any 
t-wo quantum -wavefunctions is constant in time. This dif- 
ficulty has led to the study of 'quantum chaos,' the search 
for characteristics of quantum dynamics that manifest 
themselves as chaos in the classical realm [6] [7] [8] [9] . 

As a possible signature of quantum chaos, Peres [3] [4] 
proposed comparing the evolution of a state under an 
unperturbed, H, and perturbed, H + V, Hamiltonian for 
chaotic and non-chaotic dynamics. The divergence of the 
states after a time t is measured via the overlap 

o{t) = \{Mt)\Mt))l (1) 

-where ipu is the state evolved under the unperturbed 
system operator and ipp is the state evolved under the 
perturbed operator. Recent insights have sharpened the 
differences bet-ween chaotic and regular dynamics under 
this approach and several regimes of overlap decay behav- 
ior based on perturbation strength have been identified. 
The overlap decays for a short time as a quadratic. After 
this time, for chaotic dynamics -with -weak perturbation 
the overlap decay is Gaussian, as expected from first or- 
der perturbation theory [10] [11] [12]. For stronger per- 
turbations, -where perturbation theory breaks down, the 
overlap decay is exponential. This occurs when the mag- 
nitude of a typical off diagonal element of V expressed 
in the ordered eigenbasis of H is greater than the aver- 
age level spacing of the system, A. The regime of ex- 
ponential overlap decay is called the Fermi Golden Rule 
(FGR) regime [12] [13]. The rate of the exponential de- 
cay will increase with stronger perturbation as the per- 
turbation strength squared until the decay rate reaches a 
value given by the classical Lyapunov exponent [14] [12] 
[15] or the bandwidth oi H [12]. The crossover regime 
from, Gaussian to exponential decay, has also been stud- 
ied [13]. We note that many of the works cited use 
as the fidelity. Here, we follow [11] and simply use the 
overlap, O. 

For regular, non-chaotic systems the FGR regime over- 
lap decay is a Gaussian, faster than the exponential de- 
cay of chaotic dynamics. This non-intuitive result is ex- 
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plained using a correlation functions formulation of the 
overlap by Prosen [11]. In addition, a power law decay 
oc t^/'^ has been found for an integrable system [16]. 

The initial overlap decay behavior continues until some 
saturation level [11]. For coherent and random pure 
states the saturation level 1/N for the exponential de- 
cay (in the FGR regime) and 2/N for Gaussian decay (in 
the weak perturbation regime), where N is the dimension 
of the system Hilbert space. However, for eigenstates of 
the system and mixed random states the saturation level 
increases with increasing perturbation strength [11]. 

Here we study a mixed system, a system with both 
chaotic and regular dynamics. Coherent states within the 
regular regime are practically eigenstates of the system 
and the overlap of these states oscillates close to unity 
(see figure 1) [4] [11]. This is shown in figure 1 where the 
initial coherent state is centered at a fixed point of order 
one of the regular map. Coherent states in the chaotic 
regime of the system show exponential overlap decay in 
the FGR regime and Gaussian overlap decay for weak 
perturbations. We show that in both the FGR and weak 
perturbation regimes states near the chaotic border have 
a polynomial overlap decay. 

In 1988, one of us [17] introduced in statistical mechan- 
ics the generalized entropy 
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where fc is a positive constant, pi is the probability of 
finding the system in microscopic state, i, and W is the 
number of possible microscopic states of the system; q is 
the entropic index which characterizes the degree of the 
system nonextensivity. In the limit q — > 1 we recover 
the usual Boltzmann entropy 
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Si = — fc^^Pi Inpi 



(3) 



To demonstrate that q characterizes the degree of the 
system nonextensivity it is useful to examine the Sq en- 
tropy addition rule [18]. If A and B are two independent 
systems such that the probability p{A + B) — p{A)p{B) 
the entropy of the total system Sq{A + B) is given by the 
following equation: 
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From the above equation it is realized that q < 1 corre- 
sponds to superextensivity and g > 1 to subextensivity. 
Using this entropy to generalize statistical mechanics and 
thermodynamics has helped explain many natural phe- 
nomena in a wide range of fields. 

One application of this nonextensive entropy occurs in 
one dimensional dynamical maps. As explained above, 



when the Lyapunov exponent is positive, the system dy- 
namics is strongly sensitive to initial conditions and is 
characterized as chaotic dynamics. When the Lyapunov 
exponent is zero it has been conjectured [2] (and proven 
[19] for the logistic map) that the distance between two 
initially arbitrarily close points is described by ^ 



^g.enC^"" leading to ^ = [1 -I 
{sen stands for sensitivity). 
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This requires the intro- 
duction of Ag^^„ as a generalized Lyapunov coefficient. 
The Lyapunov coefficient scales inversely with time as 
a power law instead of the characteristic exponential 
of a Lyapunov exponent. Thus, there exists a regime, 
qsen < 1, Ai =0, Xq,^^ > 0, which is weakly sensitive to 
initial conditions and is characterized by having power 
law, instead of exponential, mixing. This regime is called 
the edge of chaos. 

The polynomial overlap decay found for initial states of 
a mixed system near the chaotic border are at the 'edge of 
quantum chaos', the border between regular and chaotic 
quantum dynamics. This region is the quantum parallel 
of the region characterized classically by the generalized 
Lyapunov coefficient. 

The system studied is the quantum kicked top (QKT) 
[20] defined by the unitary operator: 



(5) 



J is the angular momentum of the top and a is the 'kick' 
strength. We use a QKT with a = 3 whose classical 
analog has a mixed phase space, regions of chaotic and 
regular dynamics. The perturbed operator used is a QKT 
with a stronger kick strength a' = 3.015. Hence, V = 
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FIG. 1. Overlap versus time for initial angular momentum 
coherent states located in the chaotic region and the regular 
region of the quantum kicked top. The system has a spin 
120 and is evolved under the kicked top Hamiltonian with 
a — Z and a' = 3.015. The overlap of the state in the chaotic 
region decreases exponentially with the number of iterations 
of the map. The state in the regular region is practically an 
eigenstate of the system and therefore oscillates close to unity. 
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The classical kicked top is a map on the unit sphere, 

+y'^ + = 1: 

x' = z 

y' = X sin{az) + y cos{az) (6) 
z' = —X cos{az) + y sin(az). 

For a = 3 there are two fixed points of order one at the 
center of the regular regions of the map. They are located 
at 

a;j = 0j = ±.6294126, y/ = .4557187. (7) 

The regular regions of the classical phase space are seen 
clearly in figure 4. 

To locate the edge of quantum chaos, we work in the oo 
(even under a 180° rotation about x and odd under a 180° 
about y] see page 359 of [4]) symmetry subspace of the 
QKT with J = 120 (in the oo subspace N = J/2). We 
set y equal to of the positive fixed point and change 
z so that the initial state can be systematically moved 
closer and further from the fixed point of the map. An 
initial state with a power law decrease of overlap is found 
at Zf — .124. The overlap decay for this state at the edge 
of quantum chaos is illustrated in figure 2 and is very 
well fit by the solution of dO/d{e) = -C^'/t^^^, (re/ 
stands for relaxation). Although we do not know how 
to derive this differential equation from first principles 
the numerical agreement is remarkable (sec also [21]). 
A time-dependent g-exponential expression analogous to 
the one shown here has recently been proved for the edge 
of chaos and other critical points of the classical logistic 
map [19]. The polynomial overlap decay is the transition 
between the quadratic and exponential overlap decays. 
This transitory region does not appear for chaotic states 
(as shown in figure 1) and is a signature of the 'edge of 
quantum chaos.' 

A power law also emerges for the above initial state in 
the weak perturbation or Gaussian regime. Here we use 
a' = 3.0003. The power law in this regime is illustrated 
in figure 2 and fit with the above equation. 
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FIG. 2. Overlap versus time for an angular momentum 
coherent state initially located at the border between regu- 
lar and chaotic zones of the QKT of spin 120 and a = 3. 
This region is called the edge of quantum chaos and shows 
the expected power law decrease in overlap. The top fig- 
ure is for a perturbation strength within the FGR regime, 
6 = .015 and the bottom figure is for a perturbation strength 
of 5 = .0003, well below the FGR regime. On the log-log 
plot the power law decay region, from about 20-70 in the 
FGR regime and 2000-7000 in the Gaussian regime, is lin- 
ear. We can fit the decrease in overlap with the expression 
[1 + {qrei - l)(i/-r9„,)^]^^'^"^'-*'^ where, in the FGR regime, 
the cntropic index Qrei ~ 5.1 and r,^^, — 40 and in the Gaus- 
sian regime qrei = 3.8 and t^^^, = 2500. The insets of both 
figures show In^^^, O = [0^~''^'' — Qrei) versus t^; since 

luq X is the inverse function of = [1 + (1 — q) x] 1-9 , this 
produces a straight line with a slope — 1/t^ (also plotted). 

The value of qrej remains constant at 3.8 for small per- 
turbations until the critical perturbation strength, 6c, 
when the typical off diagonal elements of V are larger 
than A. We can approximate 6^ — 27r /N^ = 5.4xl0~^ 
[12]. When 6 is larger than 6c, Qrei increases. The be- 
havior of Qrei and Tq^^i vcrsus 6 can be seen in figure 3. 
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FIG. 3. qrei versus perturbation strength. The value of 
qrei remains constant at 3.8 until the critical perturbation 

after which it increases. The relationship of Tq^,^^ versus per- 
turbation strength is shown on a log-log plot in the inset and 
is well fit by a line with slope -1.06. 

The location of the edge of quantum chaos for the 
QKT of spin 120 does not match up with the edge of 
the classical kicked top which appears at approximately 
Zf — .2296. This implies that classically regular regions 
of the kicked top appear chaotic on the QKT. As J is 
increased, the top becomes more and more classical and 
states exhibiting 'edge of quantum chaos' behavior are 
centered closer to the classical value for the edge of chaos. 
Hence, for J = 150, 180,210,240 the edge is observed at 
Zf - .124, .139, .151, .160 and .176, respectively. 

The edge of chaos in the quantum and classical maps 
are not obseved at the same value due to the size of the 
angular momentum coherent state. The coherent state 
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grows as J decreases causing it to 'leak out' into the 
chaotic region even though it is centered away from the 
chaotic border. This causes behavior characteristic of 
the edge of chaos to appear at different vahics depending 
on the dimension of Hilbert space. Figure 4 shows the 
wavefunctions for two values of J superimposed on the 
classical phase space. This gives an idea as to how large 
the wavefunction is compared to the regular region of the 
map. 
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FIG. 4. ClEissical phase space of the kicked top with angu- 
lar momentum coherent wavefunctions. 10000 iterations of a 
chaotic orbit starting from the point x = .6294, y = .7424, 
2 = .2294. The spherical phase space and the ellipsoidal co- 
herent states are projected onto the x — z plane (only y > 
shown) by multiplying the x and z coordinates of each point 
by R/r where R = ^/2(l - \y\) and r = - t/^) [4]. The 
regular regions of the kicked top arc clearly visible. Shown is 
a, J — 120 wavefunction (stars) and a J = 240 (circles) wave- 
function both at the edge of quantum chaos. Note that the 
variance of the J = 120 wavefunction is much larger than the 
variance of the J = 240 wavefunction. Hence, the behavior 
characteristic of the edge of quantum chaos appears further 
from the fixed point of the classical map. 

In the region of ,/ values that we have explored, no 
significant changes have been detected for qrei, because 
the changes only slightly. However, Tg^^, decreases 
with increasing J. 

To conclude, we have located a region on the border of 
chaotic and non-chaotic quantum dynamics. Quantum 
states located in this region exhibit a power-law decrease 
in overlap as opposed to the exponential overlap decay 
exhibited by fully chaotic quantum dynamics. The clas- 
sical parallel to this region is the border between regular 
and chaotic classical dynamics which is characterized by 
the generalized Lyapunov coefficient. 

The authors would like to acknowledge C. Anteneodo 
and J. Emerson for useful remarks. This work was sup- 
ported by DARPA/MTO through ARO grant DAAG55- 
97-1-0342. 



[1] A.J. Lichtonborg and M.A. Liebcrman Regular and 

Chaotic Dynamics (Springer Verlag, 1992). 
[2] C. Tsallis, A.R. Plastino, W.M. Zheng, Chaos, Solitons 

and Fractals 8, 885 (1997); M.L. Lyra, C. Tsallis, Phys. 

Rev. Lett. 80, 53 (1998); V. Latora, M. Baranger, A. 

Rapisarda, C. Tsallis, Phys. Lett. A, 273 97 (2000). 
[3] A. Peres, in Quantum Chaos, Quantum Measurement 

edited by H. A. Cerdeira, R. Ramaswamy, M. C. 

Gutzwiller, G. Casati (World Scientific, Singapore, 1991). 
[4] A. Peres, Quantum Theory: Concepts and Methods. 

Kluwer Academic Publishers (1995). 
[5] R. Schack, CM. Caves, Phys. Rev. Lett. 71, 525 (1993). 
[6] M. V. Berry, Proc. R. Soc. London Ser. A 413, 183 

(1987). 

[7] M. V. Berry in New Trends in Nuclear Collective Dynam- 
ics, edited by Y. A be, H. Horiuchi, and K Matsuyanagi 
(Springer, Berlin, 1992), p. 183. 

[8] F. Haake, Quantum Signatures of Chaos (Springer, New 
York, 1991). 

[9] O. Bohigas, M.J. Giannoni, C. Schmit, Phys. Rev. Lett. 

52, 1, 1984. 
[10] A. Peres, Phys. Rev. A 30, 1610 (1984). 
[11] T. Prosen, M. Znidaric, J. Phys. A 35 1455, 2002. 
[12] Ph. Jacquod, P.G. Silvostrov, C.W.J. Beenakker, Phys. 

Rev. E 64, 055203 (2001) 
[13] N.R. Cerruti, S. Tomsovic, Phys. Rev. Lett. 88, (2002). 
[14] R.A. Jalabert, H.M. Pastawski, Phys. Rev. Lett. 86, 2490 

(2001) ; P.M. Cucchietti, C.H. Lewenkopf, E.R. Mucciolo, 
H.M. Pastawski, R.O. Vallcjos Phys. Rev. E 65, 046209 

(2002) ; F. Cucchietti, H.M. Pastawski, D. A. Wisniacki, 
Phys. Rev. E 65, 045206 (2002). 

[15] G. Benenti, G. Casati, Phys. Rev E 65, 066205 (2002). 
[16] Ph. Jacquod, I. Adagideli, C.W.J. Beenakker quant- 
ph/0206160. 

[17] C. Tsallis, J. Stat. Phys. 52, 479 (1988). E.M.F. Cu- 
rado, C. Tsallis, J. Phys. A 24, L69 (1991) [Corrigenda: 
24, 3187 (1991) and 25, 1019 (1992)]. C. Tsallis, R.S. 
Mendes, A.R. Plastino, Physica A 261, 534 (1998). 

[18] S. Abe and Y. Okamoto, eds., Nonextenswe Statisti- 
cal Mechanics and Its Applications, Series Lecture Notes 
in Physics (Springer- Verlag, Berlin, 2001); G. Kani- 
adakis, M. Lissia, A. Rapisarda, eds., Non Extensive 
Thermodynamics and Physical Applications, Physica A 
305 (Elsevier, Amsterdam, 2002); M. Cell Mann, C. 
Tsallis, eds., Nonextensive Entropy-Interdisciplinary Ap- 
plications, (Oxford University Press, Oxford, 2002), in 
press. A complete bibliography on the subject can be 
found at http://tsallis.cat.cbpf.br/biblio.htm. 

[19] F. Baldovin, A. Robledo, PRE in press (2002); F. Bal- 
dovin, A. Robledo, Europhys. Lett, in press (2002). 

[20] F. Haake, M. Kus, R. Scharf, Z. Phys. B, 65, 381, 1987. 

[21] E.P. Borges, C. Tsallis, G.F.J. Ananos, P.M.C. Oliveira, 
cond-mat/0203348. 



4 



